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Abstract. In this paper we study the dynamics of orthogonal spatially 
homogeneous Bianchi cosmologies in R"-gravity. We construct a compact state 
space by dividing the state space into different sectors. We perform a detailed 
analysis of the cosmological behaviour in terms of the parameter n, determining all 
the equilibrium points, their stability and corresponding cosmological evolution. 
In particular, the appropriately compactified state space allows us to investigate 
static and bouncing solutions. We find no Einstein static solutions, but there do 
exist cosmologies with bounce behaviours. We also investigate the isotropisation 
of these models and find that all isotropic points are fiat Friedmann like. 



PACS numbers: 98.80.JK, 04.50.+h, 05.45.-a 

1. Introduction 

Over the past few years, there has been growing interest in higher order theories 
of gravity (HOTG). This is in part due to the fact that these theories contain extra 
curvature terms in their equations of motion, resuhing in a dynamical behaviour which 
can be different to General Relativity (GR). In particular these additional terms 
can mimic cosmological evolution which is usually associated with dark energy [1], 
dark matter [2,3] or a cosmological constant [4]. The isotropisation of anisotropic 
cosmologies can also be significantly altered by these higher-order corrections. In 
a previous work [5], the existence of an isotropic past attractor within the class of 
Bianchi type I models was found for a power law Lagrangian of the form i?". This 
feature was also found for Bianchi type I, II and IX models in quadratic theories of 
gravity [6,7]. In these cases the extra curvature terms can dominate at early times 
and consequently allow for isotropic initial conditions. This is not possible in GR, 
where the shear term dominates at early times. 

A natural extension of this analysis is to investigate the effect of spatial curvature 
on the isotropisation in HOTG. In GR, it is well-known that spatial curvature can 
source anisotropics for Bianchi models [8,9]. In this paper we extend the analysis 
in [5] to the case of orthogonal spatially homogeneous (OSH) Bianchi models [10], in 
order to investigate the effect of spatial curvature on the isotropisation of i?" models. 
OSH Bianchi models exhibit local rotational symmetry (LRS), and include the LRS 
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Bianchi types I (BI). Ill (Bill) and the Kantowski-Sachs (KS) models. For a review 
of this class of cosmologies see [10-12]). 

In GR, a cosmological constant or scalar field is required to obtain an Einstein 
static solution in a closed {k = +1) Friedmann-Lemaitrc-Robertson- Walker (FLRW) 
model [13, 14]. The existence of Godel and Einstein static universes has been 
investigated for gravitational theories derived from functions of linear and quadratic 
contractions of the Riemann curvature tensor [15]. Recently, the stability of Einstein 
static models in some /(i?)-theories of gravity was investigated [16]. It was shown 
that the modified Einstein static universe is stable under homogeneous perturbations, 
unlike its GR counterpart [18]. Static solutions are interesting in their own right, but 
are often an important first step in finding cosmologies that have a "bounce" during 
their evolution [17]. 

The existence conditions for a bounce to occur for FLRW universes in f{R)- 
gravity have been determined recently [19]. Bouncing cosmological models have been 
found for FLRW models in i?"-gravity [20,21]. This should in principle be possible for 
anisotropic models as well, since the higher order corrections can mimic a cosmological 
constant, and so prevent the model from collapsing to a singularity. In [22], it was 
shown that bounce conditions for OSH Bianchi models cannot be satisfied in GR 
with a scalar field, but can be satisfied for KS models in the Randall-Sundrum type 
Braneworld scenario. 

As in [5], we make use of the dynamical systems approach [14,23,24] in this 
analysis. This approach has been applied to study the dynamics of a range of extended 
theories of gravity [5,6,20,25-32]. However, in these works, the dynamical variables 
were non-compact, i.e. their values did not have finite bounds. This non-compactness 
of the state space has certain disadvantages (see [33] for detailed discussion of this 
issue). The standard expansion-normalised variables for example only define a 
compact state space for simple classes of ever expanding models such as the open and 
flat FLRW models and the spatially homogeneous Bianchi type I models in GR [24]. 
As soon as a wider class of models or more complicated underlying gravitational 
theories are considered, the expansion rate may pass through zero, making the state 
space non-compact (see e.g. [5,25-27]). The points at infinity then correspond to a 
vanishing Hubble factor, and the non-compact expansion-normalised state space can 
only contain the expanding (or by time-reversal collapsing) models. In order to obtain 
the full state space, one would have to carefully attempt to match the expanding and 
collapsing copies at infinity. 

While static solutions correspond to equilibrium points at infinity and can 
be analysed by performing a Poincare projection [34,35], bouncing or recoUapsing 
behaviours on the other hand are very difficult to study in this framework. In both 
cases ambiguities at infinity can easily occur, since in general only the expanding 
copy of the state space is studied. A point at infinity may for example appear as 
an attractor in the expanding non-compact analysis, even though it corresponds to a 
bounce when also including the collapsing part of the state space. 

In order to avoid these ambiguities, we will here construct compact variables that 
include both expanding and collapsing models, allowing us to study static solutions 
and bounce behaviour in i?"-theories of gravity. This approach is a generalisation 
of [13], which has been adapted to more comphcated models in [36-39]. We refer to the 
accompanying work [33] for a detailed comparison between the approach established 
here and differently constructed non-compact state spaces applied to the class of BI 
or flat FLRW models in i?"-gravity. 
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We note that we recover the isotropic past attractor found in [5] in this analysis, 
and we only obtain flat (fc ~ 0) isotropic equilibrium points. Bounce behaviour is 
found for BI, Bill and KS cosmologies, but no Einstein static solutions could be 
found in the phase space. Our analysis also reveals that we can have cosmologies that 
bounce from expansion to contraction and vice versa, depending on the value of the 
parameter n. 

The outline of this paper is as follows: In section 2 we state the field equations 
and the evolution equations for the OSH Bianchi models. In section 3 we construct a 
compact state space and then analyse the Bill and KS subspaces separately. Section 
4 is devoted to a discussion of the isotropisation of these cosmologies. 

The following conventions will be used in this paper: the metric signature is 
( — h H — h); Latin indices run from to 3; units are used in which c = SttG = 1. 



2. Preliminaries 



The general action for a /(i?)-thcory of gravity reads 

A= I dx^^fiR) + f CMdx\ (1) 



where Cm is the Lagrangian of the matter fields. The fourth order field equations can 
be obtained by varying ([1]): 

T^l = f'Rab - \fgab + Scd {g^'gab - glg'i) , (2) 

where primes denote derivatives with respect to R and Sab — VaVb/'(i?). The field 
equations ^ can be rewritten in the standard form 

Gab = Rab — \gabR = T^b'^ ' (3) 

(when f'{R) 7^ 0); where the effective stress energy momentum tensor T^l^ is given 

by 

= Kl + \gab if - fR) + - g'-^gab)] ■ (4) 

It is easy to show that the contracted Bianchi identities V^Gab = give rise to the 
conservation laws for standard matter [25] . The propagation and constraint equations 
can be obtained straightforwardly for these field equations (sec [5,40]). 

We here consider the case f{R) = R" for OSH Bianchi spacetimes, where the 
Raychaudhuri equation becomes 

e + ie^ + 2.^-i-i?-(.-i)|e + -^ = o, (5) 

and the trace free Gauss-Codazzi equation is given by 

.^-(e + (.-l)|)a+^3^. (6) 



Here O is the volume expansion which defines a length scale a along the flow lines 

„..o. e - f. .„a . ...a.a „„u„ xHe 

(see [10]). 

The Fricdmann equation is given by 



magnitude of the shear tensor is given by cr^ = ^cr^^cTat, and the 3-Ricci scalar by ^R 



le'--a'^ + in-l)^e-^-^^-^R-—^ + ^^R^O. (7) 
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Combining the Fricdmann and Raychaudhuri equations yields 

i? = 26 + fe^ + 20-2 + 3^. (8) 

We will assume standard matter to behave like a perfect fluid with barotropic index 
w, so that the conservation equation gives 

fi= -{1 + w)^lQ. (9) 

In the following, we assume n > and n ^ 1. 

3. Dynamics of OSH Bianchi cosmologies 

3.1. Construction of the compact state space 

The overall goal here is to define compact dimensionless expansion-normalised 
variables and a time variable t such that the system of propagation equations above 
©-([I]) can be converted into a system of autonomous first order differential equations. 
We choose the expansion normalised time derivative 

I d 1 d , ^ 

^Tr^Ddt 

and make the following ansatz for our set of expansion normalised variables Q: 

V3a ^^3^(i_,)^ (11) 

D ' RD^ ^ ' ' ^ 2nD^ ^ ' ' ^ ' 

3n ^ e 

K = 



Here _D is a normalisation of the form 

A, (12) 

where A is a linear combination of the terms appearing on the right hand side of the 
Friedmann equation ([7]) as discussed below. In order to maintain a monotonically 
increasing time variable, A must be chosen such that the normalisation D is real- 
valued and strictly positive. 

Note that we have chosen to define x with an opposite sign to that in [5] in 
order to have a simple form of the Friedmann equation (see below), and a can be 
both positive and negative [10]. We emphasise that the coordinates (fTTj) are strictly 
speaking only defined for R ^ 0, which means for y ^ 0. Even though the case R ~ 
may not be of physical interest, the limiting case is interesting in the context of the 
stability analysis, since we obtain equilibrium points with y = 0. This means that the 
system may evolve towards/away from that singular state if these points are attractors 
or repellers. In the analysis below we will investigate this by taking the limit y — > 
(by letting R —f 0) and find that this puts a constraint on the relation between the 
coordinates. 

We now turn to the issue of compactifying the state space. It is useful to re-write 
the Friedmann equation ([7]) as 

e^^±^-K + x + y + z = D^ + A, (13) 

where the quantities with a hat are just the variables defined in (|lip without the 
normalisation D. If all the contributions (S^, — K, x, y and z) to the central term 



X It is important to note that this choice of variables excludes GR, i.e., the case of n = 1. See [14,23] 
for the dynamical systems analysis of the corresponding cosmologies in GR. 
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in equation (|13p are non-negative, we can simply normalise with (i.e. A = 0), but 
we have to explicitly make the assumption O 7^ 0. We can then conclude that the 
state space is compact, since all the non-negative terms have to add up to 1 and are 
consequently bounded between and 1. 

However, while is always positive, —K, x, y and z may be positive or negative 
for the class of models considered herel§|. This means that the variables ([TT|) do not 
in general define a compact state space. 

In the following, we will study the class of LRS BUI models with < and 
the class of KS models with > separately, as in [13]. While we may in principle 
normalise with in the Bianchi III subspace, we have to absorb the curvature term 
into the normalisation D in the KS subspace. 

For both classes of models, we can construct a compact state space by splitting 
up the state space into different sectors according to the sign of x ,y and z. In both 
the open and the closed subspaces we will have to define 2^ = 8 sectors, corresponding 
to the possible signs of the three variables x,y, z. In the following, wc will refer to 
the spatially open Bill sectors as sector Iq to sector 80, where the subscript 'o' stands 
for 'open'. Similarly, the spatially closed KS sectors will be labeled sectors Ic - 8c, 
where the 'c' stands for 'closed'. 

After defining the appropriate normalisations for the various sectors, we derive 
the dynamical equations for the accordingly normalised variables in each sector. For 
each sector we then analyse the dynamical system in the standard way: we find the 
equilibrium points and their eigenvalues, which determine their nature for each sector. 
The overall state space is then obtained by matching the different sectors along their 
common boundaries. 

3.2. The LRS Bill subspace 

If '^R < 0, we obtain the class of spatially open LRS Bill cosmologies. This class 
of models contains the flat LRS BI models as a subclass. In this case K enters the 
Friedmann equation with a non-negative sign and does not have to be absorbed into 
the normalisation. As can be seen from the Friedmann equation in each sector (see 
Table [T]), K £ [-1, 0] and E e [-1, 1] holds in each sector. 

3.2.1. Sector Ig The first open sector denoted lo is defined to be that part of the 
state space where x ,y, z > 0. In this case all the contributions to the right-hand side 
of (|13p are non-negative, and we can choose A = 0. This means we can normalise with 
D = \Q\ = £0, where e is the sign function of O and e = ±1 for expanding/collapsing 
phases of the evolution. Note that it is crucial to include e in the normalisation: if 
we were to exclude this factor, time would decrease for the collapsing models, and 
any results about the dynamical behaviour of collapsing equilibrium points would be 
time-reversed. 

It is important to note that we have to exclude = in this sector, so we cannot 
consider static or bouncing solutions here. However, this assumption is not as strong 
as it first appears: wc can see from the Friedmann equation (jl3p that the only static 
solution in this sector appears for S: = y = z = 'S = K — 0, because all the quantities 
enter (|13p with a positive sign in this sector by construction. This means that we only 
have to exclude the static flat isotropic vacuum cosmologies!]]. Under this restriction, 

§ Note that the sign of K is preserved within the open and the closed sectors. 
II The same restriction appears in GR, see [13] 
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the normalisation above is strictly positive and thus defines a monotonically increasing 
time variable via PH)) . Equation now becomes 

1 = - A' + a; + y + z . (14) 



We can directly see from ()14p that the appropriately normalised variables (|TT|) define 
a compact subscctor of the total state space: 

X, y, ze[0,l], A'e[-1,0] and (15) 

Here Q = e is constant and not a dynamical variable. 

This sector is different from all the other sectors in both the open Bill and the 
closed KS subspaces for the following reasons. When gluing together the different 
sectors to obtain the total state space, we will actually use two copies of Iq'. one 
copy with e = 1 corresponding to expanding cosmologies and one copy with e = — 1 
corresponding to collapsing cosmologies. The two copies are in fact disconnected: The 
closed sector Ic from the KS subspace separates the expanding and collapsing copies 
of open sector Iq. Again, this refiects the fact that we cannot study static solutions in 
sector Iq. In all the other sectors we allow 9 = 0, and the expanding and collapsing 
sets are connected via the non-invariant subset (5 = 0. 

We can now derive the propagation equations for the dynamical systems variables 
in this sector by using the definitions (|10p and (jlip and substituting them into the 
original propagation equations ([S])-®. We obtain five equations, one for each of 
the dynamical variables defined in pT|) . These variables are constrained by the 
Friedmann equation (|13[) . which we use to eliminate a;, resulting in a 4-dimensional 
state space. Note that we have to verify that the constraint is propagated using all 
five (unconstrained) propagation equations, which we have done for each sector. The 
effective systerr0 is given by 



K' ^2eK 



E' = -e 

y' = \z + (2?i - 3)A' - (2n - 1)?; + (2n - + 4n - 5] 

n — 1 

n—1 ^ 



(16) 



Only in this sector does the sign of the expansion-rate appears directly in the 
dynamical equations, and we can see directly that the stability of the collapsing 
equilibrium points is given by simple time-reversal of the stability of the expanding 
points and vice versa. 

The subset K = Q (Bianchi I) is a two dimensional invariant sub-manifold, so it 
is justified to discuss the Bianchi I subspace on its own. This is done in detail in [33]. 
The vacuum subset z = and the submanifold y = are also invariant subspaces. On 
the other hand, the isotropic subset S = is not invariant unless K = 0. This agrees 
with GR, where it was found that the spatial curvature can source anisotropics for 
Bianchi models [8,9]. 

We can find the equilibrium points and the corresponding eigenvalues of the 
dynamical system (|16p . and classify the equilibrium points according to the sign of 

1 If we used the unconstrained 5— dimensional system, we would not constrain the allowed ranges of n 
and w for the different equilibrium points correctly. We would also get a fifth zero-valued eigenvalue 
for all equilibrium points. 
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their eigenvalues as attractors. repellers and saddle points (see [35]). Beeause of the 
large number of sectors that need to be studied, we do not show the results for each 
sector. Instead we combine the results from the various sectors in Table [51 

3.2.2. Sectors 2o — 8o Sectors 2o — 8o are defined according to the possible signs 
of i, z as summarised in Table [TJ In each sector A is defined as the sum of the 
strictly negative contributions to (jl3p . so that —A is strictly positive, making D 
strictly positive even for = 0. This means that 13 is a well-defined (non-zero) 
normalisation, and (jlOp defines a well-defined monotonously increasing time variable 
for each sector, even for static or bouncing solutions. With this choice of normalisation, 



Table 1. Choice of normalisation in the different LRS Bianchi III sectors, where 
the subscripts in the sector labels stand for open, differentiating the labels for 
the open sectors from the ones defined in the closed KS subspace below. We 
abbreviate x = (1 — n)RS/R, y = (I — n)R/n and z = fj,/{nR"~^). 



sector 


X 


y 


z 


normalisation 


Friedmann equation 


range of (i, y, z) 


lo 


> 


> 


> 




1 = X 4- 


-y + z + - K 


[0, 1] X [ 0, 1] 




0, 1] 


2o 


< 


> 


> 


A = i 


l = y4 


z + j:^ - K 


[-1,0] X [0, 1] 


X 


[0, 1] 


3o 


> 


< 


> 




1 = X 4 


^z + Y? -K 


[0, 1] X [-1,0] 


X 


[0, 1] 


4o 


> 


> 


< 


A^z 


1 = X 4 


^y + Y.'^ - K 


[0, 1] X [ 0, 1] 




-1,0] 


5o 


< 


< 


> 


A — X + y 


1 = z + 


Y^ - K 


[-1,0] X [-1,0] 


X 


[0, 1] 


6o 


< 


> 


< 


A^ x + z 


l = y + 


S]2 - if 


[-1,0] X [0, 1] 


X 


[-1,0] 


7o 


> 


< 


< 


A = y + z 


1 = X 4 




[0, 1] X [-1,0] 


X 


[-1,0] 


8o 


< 


< 


< 


A = X -'r y + z 


1 = S2 


- K 


[-1,0] X [-1,0] 


X 


[-1,0] 



only positive contributions remain in the Friedmann equation, and the appropriately 
normalised variables define a compact sub-sector of the total state- space, as can be 
seen from the respective versions of the Friedmann equation in Table [T] Note that 
the Friedmann equation looks different in each sector, which is of course due to the 
different normalisation for each sector. We also gain a second constraint equation 
which arises from the definition of Q: 

^ = Q'-^^ (17) 

which can be written in terms of the variables (jlip in each sector. 

It is straightforward to derive the dynamical equations for each sector, and again 
we analyse them as outlined in the previous subsection. We confirm in each sector 
that the flat LRS BI subset is indeed an invariant submanifold. 



3.2.3. Equilibrium points of the full LRS Bill state space The equilibrium points of 
the entire Bill state space are obtained by combining the equilibrium points in each 
sector. We summarise them in Table [2] Note that not all the points occur in all of 
the sectors, and some points only occur in a given sector for certain ranges of n or 
a specific equation of state w. For this reason, we cannot express all the equilibrium 
points in terms of the same variables. When possible we state the coordinates in 
terms of the dimensionless variables defined for sector lo, i.e if the given point occurs 
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in this sector. This is true for aU the points except the hne whose coordinates are 
described in terms of the variables defined in sector 2o (see below for more details on 
the relation between Ci and £2)- 

We emphasise that if the same point occurs in different sectors, it will have 
different coordinates in each of these sectors. In particular, Q can be a function of 
n or w in sectors 2o — 80 even if Q = e is a constant in sector Iq. This simply 
reflects the fact that we have to exclude the static solutions in sector lo but not in 
the other sectors. This issue will be of importance when looking for static solutions in 
section [3.4.31 In order to ensure that equilibrium points obtained in different sectors 
correspond to the same solution, we have to look at the exact solution at these points. 
This is outlined in section [3^ 

Note that each of the isolated equilibrium points has an expanding (e = 1) and a 
collapsing (e = —1) version as indicated in the labeling of the points via the subscript 
e in Table [H Similarly, the lines each have an expanding and a contracting branch 
(see below). We will however drop the subscript in the following unless we explicitly 
address an expanding or contracting solution. 

We find the three equilibrium points A, B and C corresponding to spatially flat 
Friedmann cosmologies. The expanding versions of these points correspond to the 
equally labeled points in the BI analysis [5] (see [33] for detailed comparison). These 
points were also found in the Friedmann analysis [25] . A and B are vacuum Friedmann 

points, while C represents a non- vacuum Friedmann point whose scale factor evolution 

2 

resembles the well known Friedmann-GR perfect fluid solution with a oc _ 

We now address the two lines of equilibrium points denoted by Ci and £2- Both 
these lines correspond to the spatially flat anisotropic BI cosmologies. The ratio of 
shear E and curvature component x changes as we move along both lines. We note 
that in [5] a single line of equilibrium points denoted £J was found. In section [3^ we 
will discuss in more detail how £1 and £2 are related to C\ . 

We emphasise that for Ci the two expanding and contracting branches are 
disconnected and appear as two copies Ci^^ of the line labeled by e in Tabled Each of 
these two branches range from purely shear dominated (E = 1) to isotropic (E = 0), 
to purely shear dominated with opposite orientation (E = — !)■ For £2 on the other 
hand the expanding and contracting branches are connected: each £2.+ and £2,- 
ranges from expanding (Q* > 0) and static (Q* = 0) to collapsing (Q* < 0). The two 
disconnected copies £2,+ and £2,- correspond to positive and negative values of the 
shear respectively. Note that there is no isotropic subset of £2 in analogy to the fact 
that there is no static subset of £1. 

A closer look shows that £1 and £2 are actually the same object in different 
sectors: £1 has i > hence occurs in sectors lo, 3o, 4o and 7o, while £2 is the 
analog with i < occurring in sectors 2o, 5o, 60 and 80. This statement is confirmed 
by looking at the exact solutions corresponding to the points on both lines; we find 
that both these lines have the same parametric solution of scale factor and shear (see 
section below). For this reason, we could in fact give the two lines the same label. 
However, it is useful to treat them separately, since we obtain different bifurcations 
in the sectors with a; > and a; < respectively. Furthermore, the subset of the line 
denoted by £2 allows for static solutions unlike the subset labeled £1. This is due to 
the fact that a negative curvature contribution x can effectively act as a cosmological 
constant by counter-balancing other contributions in the Friedmann equation. This 
is explored in section [3.4.31 below. 

Finally, we find the equilibrium points V and E corresponding to spatially open 
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Table 2. Equilibrium points of the full OSH Bianchi state space in terms of 
the coordinates defined for sector Iq, except for the line <C2,±, where we have to 
use the coordinates defined for Bill sector 2o (see text). Here e = ±1 labels the 

expanding/contracting solutions. We have abbreviated f{w) = ~^*'(3^ ~2)™^^^ ' 
Pi{n) = 2'n? — 271 — 1 and P2{n) = 2n^ — 5n + 5. We will not explicitly state 
the expressions for iJi (n, ui), iJ4(n, «)), which arc rational functions of n and 
w. The constants Q», S» take real values in [—1, 1]. We have denoted the 
coordinates Q that can become non-constant in sectors other than the first sector 
with the superscript f. 



Point 



{Q,K,T.,x,y,z) 



(et, 0, 0, 0) 

\ 2n— 1 ' 2n — 1 ' j 

{e\ 0, 0, iw-l, 0, 2-iw) 

/ton 3(n-l)(l+u)) (n-l)[4n-3(l+'!i))] n(13+9'!D)-2Ti^(4+3t«)-3(l+'!ij) \ 

, u, u, , ^ , ^ j 

e, 0, S,, 1-I]2^0,0) 
g., 0, ±1, Q2„ 1^0,0) 
e, -3/4, ~e/2, 0, 0, 0) 



Description 



Be 

Line £i ^ 
Line £2,± 



3u;-l 
'3iu-2' 



0, 3«;/M 



I 3(4n^-10n+7)Pi(ra) Pi(n) 6(ra''-3n+2) 9(4ra*-187i"'+3lTi^-24Tt+7) 



P2(n) 



^ Ri{n,w), eR2{n,w), ^(lLJM±!£l^ i?3(n,w), i?4(n,w)^ 



Friedmann flat 

Friedmann fiat 

Friedmann flat 

flat LRS Bianchi I 
flat LRS Bianchi I 
open LRS Bill 

open LRS Bill 
vacuum BI, Bill or KS 
BI, Bill or KS 



models. Point V is independent of n and w, while £ depends on the value of w. The 
points J- and Q can be spatially open, flat or closed depending on the value of n 
and/or w, i.e. they move through the different sectors of the total state space as n, w 
are varied. This is reflected in Tables [9] and [TOl where we summarise the stability 
properties of the equilibrium points of the closed and open subspaces separately, and 
observe that these two points occur in each subspace for certain ranges of n only. 

3.3. The Kantowski- Sachs subspace 

When "^R > 0, we obtain the class of spatially closed KS cosmologies. Here K is 
positive and needs to be absorbed into the normalisation in all sectors. This means 
that in this subspace, —A is strictly positive in all closed subsectors Ic - 8c, hence 
is strictly positive even for = 0. We can therefore consider static and bouncing 
solutions in all sectors that make up the closed subspace. The flat subspace is obtained 
in the limit "^i? — > 0. As explained in the previous subsection, we have to exclude static 
flat isotropic vacuum cosmologies in this limit. 

The closed sectors can be defined as in the Bill case, except that K no longer 
appears in the Friedmann equation (see Table [3|). Similar to the Bill case, the first 
sector labeled Ic is defined as the subset of the state space where x ,y, z > 0. In this 
case we choose A = —K{< 0), so that equation (fT3|) becomes 

l^S^+x + y + z. (18) 

The curvature can be obtained from ([TT]) . which in this sector becomes 

1^Q^ + K. (19) 
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From and it is clear that the appropriately normalised variables (fTTj) define 
a compact subscctor of the total state space with 



x,y, ze [0, 1] , K e [0, 1] and Q, S G [-1, 1] 



(20) 



Note that the variable K will not be used explicitly in any of the closed sectors. 

As in the Bill case, we derive the propagation equations for the dynamical systems 
variables in this sector and reduce the dimensionality of the state space to four by 
eliminating x via the Friedmann constraint (|18p . Again we have verified that the 
constraint is preserved using all five propagation equations. We obtain the following 
dynamical system: 



0'- o(Q -1) 



3 

^ 3Q 

y 

3g 

z 

3Q 



1 



ny 



1 



1 



2SQ(g2 _ 1 
2Sg(Q2 



ny 
n- 1 
1 



1 



+ -(S2-1)(Q2_i), (21) 



1+1 



n-1 

S2 



(l-S'-y-^) + 2g2 



ny 
n- 1 



y- z + Q'^ [1-3W + 2J:' 



—2ny 
n - 1 



We recover the following features from the Bill subspacc: The flat subset K — (here 
corresponding to = 1) is invariant, as can be seen from the Q'-equation together 
with the Friedmann equation (jlSp . Other invariant subspaces are the hyper-surfaces 
y = and z = 0. The isotropic subset E = is not invariant unless K = 0. 



Table 3. Choice of normalisation in different KS sectors, where the subscripts in 
the sector labels stand for closed. See text and caption of Table [T] for details on 
the notation used here. 



sector 


X 


y 


z 


normalisation 


Friedmann equation 


range of {x, y, z) 


Ic 


> 


> 


> 


-k 


1 = X 4 


y 4 2 4 


[0, 1] X [ 0, 1] 




0, 1] 


2c 


< 


> 


> 


A = x-k 






[-1,0] X [0, 1] 


X 


[0, 1] 


3c 


> 


< 


> 


A = y-k 


1 = a; 4 




[0, 1] X [-1,0] 


X 


[0, 1] 


4c 


> 


> 


< 


A = z-k 


1 = .X 4 


y + E^ 


[0, 1] X [ 0, 1] 


X 


-1,0] 


5c 


< 


< 


> 


A ^ X + y - k 


l = z + 




[-1,0] X [-1,0] 


X 


[0, 1] 


6c 


< 


> 


< 


A ^ X + z - k 


1 = y4 




[-1,0] X [0, 1] 


X 


[-1-0] 


7c 


> 


< 


< 


A ^ y + z - k 


1 = a; 4 




[0, 1] X [-1,0] 


X 


[-1,0] 


8c 


< 


< 


< 


A^x+y+z-k 


1 = E2 




[-1,0] X [-1,0] 


X 


[-1,0] 



The sectors 2c-8c are defined according to the possible signs of x, y, z as 
summarised in Table [3l In each sector A is defined as the sum of the strictly negative 
contributions to ^3]) . The dynamical equations analogous to l(2T|) can be derived 
straightforwardly for each sector. We then solve these equations in each sector for 
their respective equilibrium points and the corresponding eigenvalues, and classify the 
equilibrium points according to their dynamical properties. The results are combined 
with the results from the open sectors and summarised in Tables IMTUl 
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3-4-. Exact solutions corresponding to the equilibrium points 



We now derive the solutions corresponding to the various equihbrium points. Special 
attention has to be paid to the points with y = 0, since these correspond to the limit 
i? — > 0, which may makes the coordinate x singular. We will study this issue in detail 
below. Note that it is legitimate to take the limit i? ^ in the original field equations 
as long as n > 1, which results in the constraint Tj^^ 0. Consequently it is only 
possible to study the limit i? ^ for /i ^ and n > 1 when solving for the solutions 
corresponding to the equilibrium points with y = 0. 

It is important to emphasise that the dynamical system by itself is well-defined for 
y = 0; only when going back to the original equations to solve for the exact solutions 
corresponding to the equilibrium points with j/ = do we notice that there may not 
be an exact solution corresponding to these coordinates. 

We now proceed to find the exact solutions corresponding to the non-static 
{Q ,Q ^ 0) equilibrium points. As usual, we can solve the energy conservation 
equation ([9]) for the non-vacuum solutions to obtain 

^, = ^^oa-'^'+^\ (22) 

where /zq is determined by the z-coordinate of the given equilibrium point. We require 
fj,o > 0, which constrains the allowed range of n or w for a given equilibrium point 
(see below). 

In order to determine the scale factor evolution at each equilibrium point, we 
rewrite the Raychaudhuri equation (O as 

e = - (1 + g.) ^ , 



where we express the deceleration parameter qi 
dimensionless variables pT|) : 



-2^ 



(23) 

at each point in terms of the 



(24) 



Note that this equation is invariant in different sectors: for a given equilibrium point, 
each coordinate divided by is the same in all sectors. This ensures that the 
corresponding solution is invariant, no matter with which coordinates we describe 
the equilibrium point. 

Similarly, we re-write the trace free Gauss Codazzi equation ([6]) as 



1 



V3Q2 

and the curvature constraint 



El 
3 



as 



R 



1 



3. 

rt'2 



El 



(25) 



(26) 



for a given equilibrium point with coordinates [Qi 
parameter qi. 



Ki,T,i, Xi,yi, Zi) and deceleration 
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3.4-1- Power-law solutions We first study the non-stationary (g ^ —1) cosmologies, 
for which lias tlic solution 

[l + qi)t 

We have set the Big Bang time to = 0. Given 0, we can solve for all the other 
dynamical quantities for a given equilibrium point to obtain the scale factor evolution 

,-1 



ao where a = (1 + qi 



the shear 



a = — 
t 

and the curvature scalar 



const , where (3 



(1 



3- 



n2 



where 7 = 



6 



2 



1 - * + 7^ + 7^ 



(28) 



(29) 



(30) 



Again, we point out that even though a given equilibrium point formally has 
different coordinates in the different sectors, the exact solutions corresponding to the 
point are invariant, since the coordinates only enter the solutions ([28 |) -([30 | with a 
factor l/Qf. The solutions for each point are summarised in Tabled! where constants 
of integration were obtained by substituting the solutions into the original equations. 

When substituting the points with ?/ = into the original field equations, we find 
that these are only satisfied for special values of n. This is reflected in Table ID Point 
B only has a solution for n — 5/4 and w = 2/3. The solutions for points T) and 8 
only satisfy the original equations for n = 1, which has been excluded from the start. 
These points therefore do not have any physical power-law solutions. The points on 
the lines £12 only have corresponding solutions for special coordinate values, making 
only two points on each line physical (see below). 

Excluding these non-physical points, we find that the only non- vacuum solutions 
are given by C and Q. Substituting the solution ()22p into the definition of z, we find 
that the constant /io must satisfy 

2n / n \ n— 1 



Q^ 



(3n)" 



n - 1 



In order for these solutions to be physical, we require that /i > and therefore /Iq > 0. 
For C we find that this condition is satisfied for 

1 



1 < n < 



4(4 + 3w) 

while for Q it is only valid for 

1< n < , 
< 71 < iV+ , 



(13 + 9w+ V^w^ + 66u; + 73^ , w > -1, (31) 



-1 < w < 
0<w<^{- 



15 + 4V15) 



(32) 



(9 + 5w ± Vl - 'iOw - 15i 



where iV± = 

For points A and the solutions only depend on n, while the solutions at C and 
Q depend on both n and w. We can see from these solutions that points A and C are 
the isotropic analogs of points and Q respectively. 

The lines Ci and £2 have the same solutions for shear and energy density. As 
noted above, they are the same line but for different ranges of x and hence a. £1 
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contains the isotropic subset of solutions (E^ = 0) while £2 contains the static subset 
(a = Q, =0). 

In Table [5] we summarise the behaviour of the deceleration parameter q. By 
studying the deceleration parameter, we can determine whether the power law 
solutions above correspond to accelerated (— 1 < g < 0) or decelerated (g > 0) 
expansion or contraction. The expansion (or contraction) of point A is decelerating for 
n e (0, 1/2) orn € (1, i(l + V3)) and accelerating for n e (i(l + V^),2). Point B and 
lines £1,2 only admit decelerating behaviours. Point T has a decelerated behaviour for 
n£ (0, ^(1+^3)) and an accelerated behaviour for n e (i(H-\/3), 2). The equihbrium 
points C and C/ for w G [0, 1], have decelerated behaviours when n G (0, |(1 + w)) and 
accelerated behaviours when n G (|(1 + w), 00). 

3.4-2. Stationary solutions If g = —1, we obtain stationary solutions (O ~ 0), which 
have an exponentially increasing scale factor. As reflected in Table [5l the vacuum 
points A and correspond to de Sitter solutions for the bifurcation value n = 2 for 
all equations of state, while the matter points C and Q arc de Sitter-like for all n > 
but w ~ —I only, and £ appears to be de Sitter-like for w ~ I for all values of n > 0. 
Since £ has y = 0, we will have to study this case in more detail below. 
For a constant expansion rate 

e = Go , (33) 
the scale factor has the following solution 

a^aoe^^"*. (34) 
The energy conservation equation becomes 

fj, = ^ fi = fiQ . (35) 
The trace free Gauss Codazzi equation ^ can be rewritten as 



& — /3q, where /3o 



P)2 
•^0 



(36) 



3V3QI 

which on integration yields 

a^pot + ao, (37) 

where cfq is an integration constant. The evolution of the Ricci scalar can be obtained 
by substituting the solutions above into ([5]), to find 

(2+^)e2 + 2(/3ot + ao)^ (38) 

As before, we substitute the solutions at each equilibrium point into the definition 
of the coordinates, which constrains the constants of integration for each point. In 
particular, /3o = holds for all stationary equilibrium points, which means that we 
only have constant or vanishing shear. 

As in the power-law case, we see that all the equilibrium points except for C and 
Q correspond to vacuum solutions fi ^ 0. For point C the energy density is given by 

^^^^c ^4"-i3-"(2-n)e^", (39) 

and for Q the energy density is given by 

/i = M?=4"-^(2-n)eg". (40) 
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Both of these solutions only hold for 1 < n < 2 with w = — 1. 

Again, we substitute the generic solutions into the original field equations for each 
point, and find that the original equations are satisfied for all points with y 0- It is 
however not possible to find a stationary solution at point £ (which has y = 0), even 
after carefully considering the limit y — > 0. 

3.4-3. Static solutions The static equilibrium points arc characterised by = O = 0. 
These points satisfy Q = x = 0, where the second identity comes from the fact that 
if Q = 0, then we require that a; = from the definition of the variables, as discussed 
below. 

We will now explore which of the equilibrium points obtained above correspond 
to static solutions. As indicated above, even though Q = ±e holds in the first sector 
as stated in Table [21 Q can be a function of n and/or w in the other sectors. In order 
to find the static equilibrium points, we have to look at the coordinates that each 
equilibrium point takes in each sector, and find the values of n and/or w for which 
Q = in the given sector. 

An obvious static solution appears to be the subset Q = on line C2,± for all 
values of n and w. We can however not find a solution corresponding to this limit, 
since Q* = implies cr = 0, which contradicts the value of the shear coordinate of 
this equilibrium point. We can study the eigenvalues associated with the line C2 in 
the limit Q — *■ 0^ and find that the static subset is an unstable saddle point for all 
values of n for both £2,+ and £2,-- 

The point A appears to admit a static solution for the bifurcation value n = 1/2. 
This bifurcation only occurs in sectors 2, 3, 6 and 7 of the open and the closed sectors. 
However, it is not possible to find a solution satisfying the coordinates of the static 
equilibrium point that satisfies the original field equations. For this reason, this static 
equilibrium point is unphysical. We explore the stability of the static solution in the 
limit n — > 1/2 from the appropriate sides: for example, point A only lies in the open 
sector 2 for n e [0, 1/2] or n e [2, 00], making only the limit n 1/2" well-defined. 
We find that this bifurcation represents a saddle point in the state space since two of 
the eigenvalues approach 00 from the left and —00 from the right, making the point 
unstable. 

Even though the Q-coordinate of point S is a function of w in sectors 2, 4, 5 and 
7, Q cannot be zero for any values of w. This means point B does not admit any static 
solutions. 

Point C can only be static in the limit n — > in sector 6forw = 0, 1/3, 1 and in 
sector 3 and 5 for w = — 1. Again, we cannot find a solution for this special case, but 
this case is physically not interesting either way. 

The (5~coordinate of point £ is zero in sectors 6-8 for w = 2/3, but again there 
is no solution corresponding to this limit. 

Even though point J- has Q as a function of n in open sectors 2 and 6, Q{n) is 
non-zero for the allowed ranges of n. 

Point Q becomes static in the limit n —^ in sectors 4, 6 and 8, which again is 
not physically relevant. 

* Note that unlike in the bouncing or recoUapsing case below, we do not consider Q = y = 0, x ^ 
here, since this corresponds to the limit i? — > 0. While we may want to study a bounce where the 
Ricci scalar approaches zero and then grows again, we arc not interested in static solutions that have 
vanishing Ricci curvature at all times. 
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Table 4. Solutions for scale factor, shear, curvature and energy density 
corresponding to the equilibrium points. 



Point 



Scale factor (a) Shear (ct) 



Ricci Scalar (i?) 



A 
B 
C 

T 
Q 

Line 
Cx 



w 



n^2 
n = 2 

_ 5 
4 ' 

-1 

w = —1 
n = 2 



w ^ 
w = 



■^2 _ 5-4n 



2ri-l ' 

2k -1 
5— 4n 



«G (1,1) 

, n e (i,i) 



aQ\t\~ 

ao|i|^ 
ao|i|3(i+™) 



ao\t\' 
ao|t|3(i+«) 



(l+2ra-2ra^) 
V3(n-2)t 

(3(l+u;)-27t) 
V3(l+iu)t 

V3|S»| 
(2+S^)t 

x/3|Q,| 
(l+2QJ)t 



6n(l-n)(2n-l)(4rt-5) 
(ri-2)2t2 



4K(4ra-3(l+'U))) 

3(l+u>)2t2 
4C)2 
3*^0 

6n(l-ra)(4ra''-10rt+7) 
{n-2)^t^ 

10o2 
3 

4n(2(ra-l)+-i«(l+3m-2rt)) 
(l+«))2t2 

3 *-^0 








Table 5. Deceleration parameter for the equilibrium points. In the last three 
columns we state explicitly for which values of n the deceleration parameter q 
(stated in the second column) is less, equal to or larger than 0, i.e. whether the 
have accelerated, de— Sitter-like or decelerated behaviours. The parameters are 



P+ = 1(1 + ^3) and Sv_ 



\{l + w). 



Point 



q=-l -1< g < 



Range of n 



g > 



A 


l+2n-2n^ 
l-3n+2ra^ 


all 


2 




/ (0,1) 

I (1,^^+) 


B 


1 


all 






(0,(X3) 


C 


3(l+io)-2n 


-1 


(0,00) 






In 


[0,1] 




(S'tu,oo) 


(0, S^) 


T 


l+2n-2n^ 
5-5n+2n^ 


all 


(0,00) 


(^+,2) 


(0,P+) 


Q 


3(l+iD)-2n 


-1 








2n 


[0,1] 




(5't„,oo) 


(0, S^) 


Line 












£1 


1 + 


all 




(0,cx)) 




£2 


1 + ^ 


all 




(0,00) 
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5.5. The full state space 

The full state space is obtained by matching the various sectors along their common 
boundaries. Because the full state space is 4-dimensional it is not easily visualised, 
so we refer to [33] for an illustration of the 2-dimensional Bianchi I vacuum subspace 
and the 2-dimensional flat FLRW subspace with matter. We emphasise that we have 
to formally exclude the subset with Q = Q and x ^ Q from the state space unless 
y ^ 0. This is an artifact of the definition of the variable cc, and reflected by the fact 
that there are no orbits crossing this subset - the only trajectories crossing the plane 
Q ~ Q pass through the points with a; = or j/ = 0. 

3.6. Qualitative Analysis 

We summarise the dynamical behaviour of the equilibrium points and lines of 
equilibrium points in Tables [SI [S and [TO] and Tables [7| and [5| respectively. For the 
stability analysis, we only consider the four cases: cosmological constant w = — 1, dust 
w = radiation u; = 1/3 and stiff matter w = 1. We only state the results for the 
equilibrium points corresponding to expanding solutions. The collapsing points are 
obtained by time-reversal - in other words their dynamical stability properties arc 
simply reversed: If A+ is a rcpellcr for a given range of n, then A- is an attractor for 
the same range of n. 

Table [6] consists of all the BI subspace equilibrium points (excluding the lines) ; 
their behaviour is similar to the flat Friedmann points which were found previously 
[5,25]. We note that some of the solutions corresponding to the Friedmann and BI 
equilibrium points, have been found in [41,42]. 

The lines of equilibrium points have to be treated more carefully Wc summarise 
their dynamical behaviour in Tables [7| and [51 As noted above, the two lines include 
the same parametric solutions, but Ci corresponds to x > 0, while £2 has x < 0. 

Since these lines correspond to flat solutions they should have been found in [5]. 
In fact, the authors of [5] found a line of equilibrium points denoted by extending 
over S** S [0, 00), where S** measures the shear contribution to the Friedmann 
equation. The range < S** < 1 corresponds to our £1.+ , while T,^^ > 1 corresponds 
to £2,+ as can be seen from the solutions of the scale factor in [5]. In [5] the isotropic 
solution was at S** = and the static one occured for S** — > 00. Note that [5] did 
not address the collapsing solutions £1,- or opposite orientation of the shear £2.- 
because the phase space for BI is symmetric about the plane E = 0. The stability for 
I] < can be obtained by time reversal from the corresponding points in the S > 
subspace. 

As stated above, £1.2 arc only physical for certain special coordinates. These are 

/ 5 — An /2n — 1 

S*-±y^^^, ne (1,5/4), and =±^^-^, n e (1/2, 1), 

for £1 and £2 respectively. These equilibrium points are always saddles in nature. 

The nature of the Bill equilibrium points is stated in Table [9l For the sake of 
completeness we have included the stability of points V and £, but will not discuss 
them any further since they are not physical. The point lies in the Bill subspace 
for n g (0,i(l + V3)). is a saddle for u; = — 1 and for n G (1, 5/4) when w = 0, 
but an attractor otherwise. G lies in the Bill subspace for n E (0,3/2) when w = 0, 
for all n when w = 1/3 and n € (0, 1) and n € (3, 00) when w = 1. G+ is saddle except 
for n G (1, 5/4) when w = where it becomes an attractor. 
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Table 6. Nature of the expanding (e = +1) spatially flat BI equilibrium 
points. The collapsing analogs arc simply time reversed. The parameters are 
P+ = 1(1 + V3) and V+ = jij(ll + v^). 



Point 


w 








Ran^ 


TQ of n 












fO 1/2) 


fl/2 1) 


fl V-l) 


(Vs. 5/4~) 


('5/4 P^) 


(P^ 3/2) 


(3/2 2) 


(2, oo) 


A+ 


-1 


Saddle 


Attractor 


Repeller 


Repeller 


Saddle 


Saddle 


Saddle 


Attractor 







Saddle 


Attractor 


Repeller 


Repeller 


Saddle 


Attractor 


Attractor 


Attractor 




1/3 


Saddle 


Attractor 


Repeller 


Repeller 


Saddle 


Attractor 


Attractor 


Attractor 




1 


Saddle 


Attractor 


Repeller 


Saddle 


Saddle 


Attractor 


Attractor 


Attractor 


B+ 


-1 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 







Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 




1/3 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 




1 


Repeller 


Repeller 


Saddle 


Saddle 


Saddle 


Saddle 


Repeller 


Repeller 


c+ 


-1 


Saddle 


Saddle 


Attractor 


Attractor 


Attractor 


Attractor 


Attractor 


Saddle 







Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 




1/3 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 




1 


Saddle 


Saddle 


Saddle 


Repeller 


Repeller 


Repeller 


Saddle 


Saddle 



Table 7. Nature of the line of expanding spatially fiat anisotropic equilibrium 
points Here Si,(n) = y 2;T-l ^ bifurcation value depending on n. 



w 


n 


S e hl,-Sb(n)) 


S e (-Eb(n),E6(n)) 


S e {Sb(n),l) 




ne (0,1) 


Repeller 


Repeller 


Repeller 


-1,0, 1/3 


ne (1,5/4) 


Repeller 


Saddle 


Repeller 




n > 5/4 


Repeller 


Repeller 


Repeller 


1 


All n 


Saddle 


Saddle 


Saddle 



The nature of the KS equiUbrium points is stated in Table [TOl Point hes in the 
KS subspace for n > ^(1 + -y/S) and !F+ is always a saddle. Similarly, Q lies in the KS 
subspace for all n when w = — 1, for n > 3/2 when w ~ and n G (1,3) when w = 1. 
is saddle except for n £ (1, 1.13) when w = 1, where it is a repeller. 
We can identify the following global attractors and repellers: A+ is a global 
attractor for G (P+,2) when w ~ 0, 1/3 and 1, and for n G (2, oo) (all w). When 
w = —1, C+ is a global attractor for n e (1, 2) and for n S (0, 1/2). Point !F+ is 
a global attractor for n G (0, 1/2) and n E (5/4, P+) when w = 0, 1/3 and 1, and for 
n G (1,5/4) when w = 1/3 and 1. t/+ is only a global attractor for 7i e (1, 5/4) when 
w = 0. By time reversal the corresponding contracting solutions are global repellers. 
There are no global repellers in the expanding subspace since the lines Ci^2 contain 
repellers and hence there are no global attractors in the collapsing subspace. 

3. 7. Bouncing or recollapsing trajectories 

As motivated above, any trajectory corresponding to a bouncing or recollapsing 
solution must pass through x=^Q^Qovy = Q = Q. 
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Table 8. Nature of the line of spatially flat anisotropic equilibrium points 
C2,±- Here Qh(n) = ^ 5" 4^ is ^ bifurcation value depending on n. We discuss 
the bifurcation Q = in the section on static solutions below. Note that the 
dynamical behaviour of C2,+ and C2,- is identical. 



n 


Q e l-l,-Qt{n)) 


Q e (-Qb(n),o) 


Q e (0,Q6(n)) 


Q e (Qb(n),i] 


n 6 [0,f/2] 


Attractor 


Attractor 


Repeller 


Repeller 


n e (1/2, 1) 


Attractor 


Saddle 


Saddle 


Repeller 


n > 1 


Attractor 


Attractor 


Repeller 


Repeller 



Table 9. Nature of the spatially open Bianchi III equilibrium points, where 
P+ = 1(1 + V3). 



Point 


w 






range of 


n 










(0, 1) 


(1,5/4) 


(5/4, P+) 


(^^+,3/2) 


(3/2,3) 


(3,00) 




All 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


e+ 


-1 


Attractor 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 




0, 1/3, 1 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


^+ 


-1 


Saddle 


Saddle 


Saddle 













Attractor 


Saddle 


Attractor 










1/3 


Attractor 


Attractor 


Attractor 










1 


Attractor 


Attractor 


Attractor 








Q+ 


-1 



















Saddle 


Attractor 


Saddle 


Saddle 








1/3 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 


Saddle 




1 


Saddle 










Saddle 



The existence of bouncing orbits for Bianchi I models has been studied in [33]. 
In the vacuum case it was found that there exist bouncing/recoUapsing trajectories, 
but only for y < 0. If 71 > 1, i? has to be negative and there can only be re-collapse 
(G) < 0). For n g [0, 1/2] re-collapse may occur if i? > 0, and for n G [0, 1/2] there may 
be a bounce {Q > 0) for positive R. In all cases, the bouncing trajectories have to pass 
through the single point x = Q = (denoted by Ai in [33]) in the 2-dimensional BI 
vacuum subspace. Note that it is not possible to achieve a bounce through y = Q = 
here, since a line of equilibrium points passes through that point in this subspace. 

When matter is added, we obtain another degree of freedom, and unlike in GR, 
the matter term may enhance bouncing or recoUapsing behaviour due to the 
term coupled to the energy density. The corresponding trajectories now have to pass 
through the 1-dimensional lines with x = Q = Q oy y = Q = instead of the single 
point M. 

In the presence of spatial curvature, it is yet easier to achieve bouncing or 
recollapsing behaviour. If < 0, the results from the flat Bianchi I case are 
qualitatively recovered. For > however, there are differences to the Bianchi 
I case. In particular, positive spatial curvature allows = even for positive y. 
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Table 10. Nature of the spatially closed Kantowski-Sachs equilibrium points, 
where P+ = i(l + \/3) and X 1.13. 



Point 


w 




Ranf 


;e of n 










(0,1) (1,X) 




(f+,3/2) 


(3/2,3) 


(3, oo) 


^+ 


-1 






Saddle 


Saddle 


Saddle 











Saddle 


Saddle 


Saddle 




1/3 






Saddle 


Saddle 


Saddle 




1 






Saddle 
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4. Isotropisation of OSH Bianchi models 

It is possible to study isotropisation by looking at the stability of the Friedmann 
points in the state space (see [24] and references therein). When such an isotropic 
point is an attractor, then we have asymptotic isotropisation in the future. If the 
point is a repeller we have an isotropic initial singularity, and when it is a saddle 
we have intermediate isotropisation. Because of the dimensionality and complexity 
of the state space, we will not study specific orbits to investigate viable models. We 
will therefore restrict the following discussion to the behaviour around the equilibrium 
points only. 

In the previous sections we found two isotropic equilibrium points that admit 
cosmological solutions: a vacuum point A and a non- vacuum point C. These points 
were also found in the BI case [5]. 

The expanding point A+ is an isotropic past attractor for n G (1,5/4) when 
w = —1, or 1/3, and for n £ (1, jij;(ll + a/37)) when w ^ 1. As pointed out 
in [5], this is an interesting feature, since the existence of an isotropic past attractor 
implies that we do not require special initial conditions for inflation to take place. The 
contracting analog A- is an isotropic future attractor in these ranges. A+ is a future 
attractor for n > 2 when w = —1 and for n > ^(1 + VS) when u; = 0, 1/3 or 1. By 
time reversal. A- is a past attractor for these parameter values. 

The equilibrium point C+ is an isotropic past attractor for n G ( ^-^(li -|- ^/37, 3/2) 
when w = I and an isotropic future attractor for n € (1,2) when w = —1. When 
w = or w = 1/3 this point is a saddle for all values of n. This means that in this 
case we have a transient matter/ radiation dominated phase in which the model is 
highly isotropic and hence potentially compatible with observations. 

We note that all isotropic equilibrium points found in this analysis are flat 
Friedmann like, unlike in [25], where the isotropic points A and C with non-zero 
spatial curvature were found. The reason for this discrepancy is that the plane E = 
is no longer invariant when allowing for non-zero spatial curvature (fc ^ 0); as in 
GR spatial curvature causes anisotropies to grow in models with i?"-gravity. For this 
reason the points A and C no longer remain equilibrium points in the full OSH Bianchi 
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state space. 

There are two equilibrium points of interest with non-zero shear: the vacuum 
point J- and the non- vacuum point Q. These points are isotropic for certain bifurcation 
values of n and w: T \s isotropic for n = ^(1 + V?>) for all w, and Q is isotropic for 
n = |(1 + w) \i w > —1. The KS point is a past attractor for n E (1,1.13) 
when w = 1 and a saddle for n > 3/2 when w = 0. This means that we can 
have initial conditions which are anisotropic, or we can have intermediate anisotropic 
conditions which are conducive for structure formation, provided that the anisotropics 
are sufficiently small. When w — 1/3, the point lies in the Bill state space and 
is a saddle for all values of n, and when w = the same applies for n € (0, 1) or 
n e (5/4,3/2). 

5. Remarks and Conclusions 

Our main aim in this paper was to investigate the effects of spatial curvature on the 
isotropisation of OSH Bianchi models in i?"-gravity, and to possibly identify static 
solutions and bounce behaviours. To achieve this goal, we constructed a compact state 
space which allows one to obtain a complete picture of the cosmological behaviour for 
expanding, contracting and static as well as bouncing or rccollapsing models. This is 
not possible with the non-compact variables used in [5], since the equilibrium points 
with static solutions do not have finite coordinates in this framework. The Poincare 
projection also does not allow one to patch together the expanding and contracting 
copies of the state space, so bounce behaviour cannot be investigated. This is discussed 
in detail in [33] for the BI subspace, where the results obtained here are compared to 
the results obtained in [5]. 

We do not find any exact Einstein static solutions in this analysis. However 
we do find orbits that exhibit cyclic behaviour, which was expected from previous 
work examining the conditions for bouncing solutions in f{R) gravity [19]. We also 
recover all the isotropic equilibrium points that were found in [5]. The expanding 
vacuum point A+ is a past attractor for n S (1, 5/4) as in the BI case. We emphasise 
that we only find flat (k = 0) isotropic equilibrium points {A, B and C). Therefore 
for these types of theories, isotropisation also implies cosmological behaviours which 
evolve towards spatially flat spacetimes. Late time behaviour with non-zero spatial 
curvature will have a growth in anisotropics, as in GR. 

In conclusion, we have shown that spatial curvature docs indeed affect the 
isotropisation of cosmological models in i?"-gravity. While no exact static solutions 
could be found, we did find that bounces can occur in these cosmologies. 
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